An Eshelby model for the highly viscous flow 
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The shear flow and the dielectric a-process in molecular glass formers is modeled in terms of 
local structural rearrangements which reverse a strong local shear. Using Eshelby's solution of the 
corresponding elasticity theory problem (J. D. Eshelby, Proc. Roy. Soc. A241, 376 (1957)), one can 
calculate the recoverable compliance and estimate the lifetime of the symmetric double- well potential 
characterizing such a structural rearrangement. A full modeling of the shear relaxation spectra 
requires an additional parametrization of the barrier density of these structural rearrangements. 
The dielectric relaxation spectrum can be described as a folding of these relaxations with the Debye 
process. 

PACS numbers: 64.70.Pf, 77.22. Gm 



I. INTRODUCTION 

From an experimental point of view, broadband dielec- 
tric spectroscopy-^ is the most versatile method to study 
the flow process in undercooled molecular liquids. How- 
ever, the relation between dielectric relaxation and shear 
flow is not yet clear. The classical Debye picture and its 
extension to viscoelasticity^ considers the molecule as a 
small sphere immersed in the viscoelastic liquid. It pre- 
dicts a slow dielectric decay, about a factor of ten slower 
than the one found in experiment^. A thorough quanti- 
tative analysis of dielectric and shear data in seven glass 
formers® showed a general qualitative agreement with the 
extended Debye scheme, but a rather poor quantitative 
fit. 

One cannot help feeling that the extended Debye 
scheme mistreats the structural rearrangements of the 
highly viscous fluid. The Debye relaxation time of the 
molecular orientation is usually longer than the Maxwell 
relaxation time of the shear stress, the more so the larger 
the molecular volume is. No such retardation is expected 
for a local structural rearrangement, which ought to be 
characterized by the same relaxation time for shear and 
dielectrics. Thus the modeling should rather separate 
the viscous effects from those of the structural rearrange- 
ments, both in shear and dielectrics. Also, in the present 
unsettled state of understanding of the highly viscous 
flow, with many different ideas and recipes^&^ i 1 ^ 11 , an 
attempt to understand the viscous flow itself in terms of 
a sequence of structural rearrangements in time seems 
legitimate. This is the purpose of the present paper. 

In order to contribute effectively to the flow, the struc- 
tural rearrangement should change the shape of the re- 
arranged region in the direction of the flow. This implies 
a strained state of the embedding matrix, against the 
flow direction before the jump and in flow direction after. 
This mechanism will be explained in detail in the next 
section, section II. It leads to a finite lifetime of the cor- 
responding double-well potential, because the surround- 
ing matrix is itself able to flow. This lifetime implies 
a specific cutoff function for the barrier density of the 



structural rearrangements. The comparison with experi- 
ment in section III requires a specification of the barrier 
density, which does not follow immediately from the pic- 
ture. As we will see, this requires three parameters even 
if there is no Johari-Goldstein secondary relaxation peak, 
and three more if such a peak is present. Such a large 
number of parameters makes it difficult to check the va- 
lidity of the model with any certainty from the rather 
broad and featureless relaxation spectra. 



II. ESHELBY MODEL FOR THE HIGHLY 
VISCOUS FLOW 

A. Shear strain defects 

The central concept of this model is a structural re- 
arrangement of a limited region in the sample which 
changes its shape to a sheared one. A simple example 
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FIG. 1: Elementary flow process (schematic), involving a rear- 
rangement of four closely packed spherical atoms or molecules. 
The dashed lines show the shape of the hole containing the 
four molecules in the unstrained matrix, the continuous lines 
the shape of the hole in the strained matrix. 
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of such a rearrangement is the special case of four closely 
packed spherical atoms or molecules shown in Fig. 1. 

Imagine that the flow occurs by expansion in the verti- 
cal direction and a contraction of the horizontal direction 
of the figure. After a while, the four molecules at the top 
of the figure will no longer experience an adapted ma- 
trix, but rather a matrix which tends to compress them 
in the horizontal direction and to tear them apart in the 
vertical one. Since the four molecules have the alterna- 
tive of the equivalent close packing shown at the bottom 
of the figure, there is a moment in time at which one 
has a symmetric double-well potential for molecules and 
matrix, the situation of the center of the figure. 

The Eshelby model described here supposes that the 
elementary flow process in molecular liquids is of this na- 
ture, not necessarily restricted to four molecules (in fact, 
the quantitative discussion of section IV will show that 
the elementary flow process must be more complicated) , 
but involving a structural rearrangement by a thermally 
activated jump between two stable locally ordered con- 
figurations with opposite shear strain with respect to the 
embedding matrix. 

The physical problem of a small piece of matter able to 
transform to a sheared shape within an elastic matrix has 
been treated fifty years ago by several authors, notably 
by J. D. EshelbyiSiis. 

Here, we translate Eshelby's result into the usual con- 
vention, in which the shear angle e and the shear stress a 
are related by a — Ge (G infinite frequency shear modu- 
lus) at short times and by er = rje (rj viscosity) in the long 
time limit. Let v be the volume of the spherical inclusion 
and ei the shear angle difference between its two stable 
configurations (in the example of Fig. 1, = 7r/3). Then 
the energy of the two equally strained configurations of 
the center of Fig. 1 is 
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The coefficient 7 is given by 



1 = 



7 — 5erp 
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where ap is Poisson's ratio. Since Poisson's ratio lies 
between 0.1 and 0.33 for the known glasses, 7 lies between 
0.48 and 0.533, close to 1/2. 

Eshelby's solution divides the energy into two almost 
equal parts, one located in the inclusion and one outside. 
Their ratio is 7/(1 — 7). The inclusion would have to 
distort by e^/2 to fit exactly into the unstrained hole. In 
the four-atom case of Fig. 1, this is the saddle point en- 
ergy for the thermally activated jump between the two 
stable configurations. It is considerably less than twice 
the energy of the two stable Eshelby minima because the 
saddle point has a lower energy than the harmonic ex- 
trapolation. Thus, in our simple four-atom example, the 
barrier height V is lower than the formation energy E a 
of the symmetric double-well potential. We will see later 



that this cannot be true for the defects which actually 
destabilize the amorphous solid. 

The Eshelby treatment supplies as well the interaction 
energy of the strain defect with the external stress com- 
ponent a oriented along the shear strain 



avei 



(3) 



This implies that the asymmetry A between the two min- 
ima changes by avei/2 in the presence of an external 
stress a, or by eGvei/2 in the presence of an applied 
external strain e. 



B. Double- well potential lifetime and cutoff 
function 



The lifetime r c of such a strain defect in the viscous 
liquid is estimated in a simple mean-field consideration. 
The shear stress outside, where other and independent 
strain defects exist, can decay according to the Maxwell 
relaxation time 
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We assume that the energy inside the inclusion can only 
decay by the yielding of the surrounding matrix. Then 
the energy which has to decay is about twice the outside 
energy, which implies that r c is about twice r m . More 
accurately, one has to take into account that the outside 
energy is not pure shear; about ten percent of it is com- 
pressional energy, which cannot be expected to decay. 
Thus one should have 



(5) 



The double-well potential lifetime r c determines the 
cutoff-barrier V c of the strain defects participating in the 
flow process according to the Arrhenius equation 



tv = r exp(]//fc B T), 



(6) 



where T is the temperature and tq is a microscopic life- 
time of 10~ 13 s. 

At the critical energy barrier V c , the potential decay 
begins to become faster than the jumps over the barrier. 
If the barrier is higher than V c , the double- well poten- 
tial begins to flow away before the population of the two 
minima can equilibrate. Therefore these higher barriers 
will not participate in the flow process. 

Since one has two competing equilibration processes, 
the potential decay with the time constant t c and the 
thermally activated jumps over the barrier with the time 
constant ry, the contribution of the jump mechanism is 
given by the cutoff function l c (V) 



1 



t c + t v 1 + cxp((^ - V c )/k B T) 



(7) 
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C. Stationary flow and recoverable compliance 



III. COMPARISON TO EXPERIMENT 



In a stationary flow e, the double-well potential asym- 
metry A changes continuously. If we ascribe e to the 
continuous passage of inclusions from one stable config- 
uration to the other, it must result from an integral over 
all these processes. The constant flow e can be deter- 
mined by counting all strain defects passing through the 
asymmetry zero 



dV 



G{vei) 2 n(V,0) 



lc(V), 



(8) 



The factor 1/5 results from the averaging over the dif- 
ferent strain orientations and l c (V) is the cutoff at the 
barrier V c from eq. ([7]). Thus 



dV 



(ve^ 2 n(V,0) 
20 



(9) 



Equation ([9]) can be used to calculate the contribution 
of the shear strain defects to the shear compliance. To 
do this, consider first the free energy F — —ksTlnZ of 
a single strain defect with asymmetry A. The partition 
function Z is 



2 cosh ■ 



A 



2k B T 



(10) 



Since the asymmetry A changes by avei/2 if one applies 
a stress a in the direction of e<, one has a contribution 
to the shear compliance determined by 



d 2 F 



Wk B T cosYf A/2k B T 



(11) 



In order to get the full recoverable compliance J® 
(the compliance after subtraction of the viscous flow 
contribution 14 ), one has to integrate over the barrier 
heights V and the asymmetries A. Here, we assume that 
the dependence of n(V, A) on A is given by the Boltz- 
mann factor exp(— F/k B T) 



n(V, A) = n(V, 0) cosh - 



A 



2k B r 

The equation for the recoverable compliance is 
1 r [°° v 2 e 2 n(V,0) 
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dVdA 



80k B TcoshA/2k B T' 



(12) 



(13) 



The asymmetry can be integrated out. With equation 
(J5J, one finds for the relaxational part of the recoverable 
compliance 



GJ° e -l = f 



(14) 



A. The parametrization of the barrier density 

The comparison to experiment requires the knowledge 
of the barrier density n(V, 0), which is difficult to model, 
particularly for glassforming systems, which need to be 
complex to avoid crystallization. Therefore we describe 
n(V,0) or, more accurately, its product with the average 
square of the coupling constant, by the pragmatic form 



l(V) 
Uc(V) 



exp ■ 



P(V - V c ) 



(V - v c ) 



k B T 



1 

— exp ■ 

24 v 5k B T 



(15) 



The normalization factor l r . 
ization condition 



is determined by the normal- 



dVl(V) = 1. 



(16) 



The parameter (3 has a close correspondence to the 
Kohlrausch-/? of approximately 0.5 of the Kohlrausch 
form exp — 1@ , the most popular fitting form for glassy re- 
laxation. It supplies the slope lo~^ at frequencies shortly 
above the a-peak frequency. The second term in eq. (|15[) 
provides a small curvature of this slope toward very high 
frequencies; we will see that one needs it. The prefac- 
tor of this term and the exponent are in fact additional 
parameters, which were fitted to the dielectric data of 
propylene carbonate 1 ^ in Fig. 2. But the form is purely 
pragmatic, without any physical significance, with the 
exception of the cutoff. Note that this cutoff is not 
very sharp; therefore the /3-values one obtains tend to 
be larger than those of a Kohlrausch fit, between 0.5 and 
0.7. 

If the glass former exhibits a pronounced secondary 
peak (Johari-Goldstein peak) one needs to add a gaussian 
with three parameters to the form of eq. (fT5]) . In this 
case, one has to reckon with different peak amplitudes in 
shear and dielectrics 1 ^. 
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The comparison to experiment in section III shows larger 
values. This will be discussed in section IV. 



FIG. 2: The dielectric data^ of propylene carbonate at 158 
K to which the barrier density was adapted. 
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The complex shear compliance, from which the com- 
plex modulus G(lo) can be easily calculated by inversion, 
is given by 



GJ'H = 1 + (GJ e ° - 1) / l(V) 



1 



and 



GJ» = (GJ e °-l) / l(V) 



LOTy 



1 



-dV- 



dV (17) 



(18) 



Note that here appears the factor fo = GJ® — 1, which 
according to eq. (|14[) should be 7r/2. r m is the Maxwell 
time r\jG. 

If one deals with a single type of double- well potentials, 
they should show up with the same l(V) in the shear 
compliance and in the dielectric susceptibility. The only 
difference lies in the viscous effects, which in the dielectric 
case should lead to a relaxation with the Debye relaxation 
time td 



td 



k B T • 



(19) 



where rn is the hydrodynamic radius of the molecule. 

Since structural rearrangements and Debye relaxation 
must be considered as independent processes influencing 
the same molecule, one has to fold the two processes in 
frequency or to multiply them in time. Thus one can use 
the same l(V) as in the shear case, replacing the ry of 
the Arrhenius equation, eq. ([5]), by the shorter relaxation 
time r„ 



TyT D 

T V +T D ' 



(20) 



With this definition, the dielectric susceptibility (after 
subtraction of the conductivity contribution) reads 



and 



e'(uj) - e c 
e(0) - eoc 



e(0) - e c 



l(V) 



1 



1 + W 2 T,2 



dV 



l(V)- ^dV. 



(21) 



(22) 



Here e(0) is the static dielectric susceptibility, is the 
real part of e(u>) in the GHz range (larger than n 2 , the 
square of the refractive index, because of vibrational 
contributions^) . 

These equations for the real and imaginary part of the 
dielectric constant contain one implicit assumption which 
might be wrong: That the same double-well potentials 
which allow the sample to flow are able to fully equili- 
brate the molecular orientation. We will come back to 
this point in the discussion. 



B. Examples 

The first example is Santovac 5P or PPE (polypheny- 
lene), a diffusion pump oil, consisting of a short chain of 
five phenyl rings connected by oxygens. The measure- 
ments were supplied by the Roskilde group£. 

If one adjusts the five parameters G, T], J®, /3 and r c to 
the shear data, they exhibit a strong scatter; the data are 
too featureless to be able to fix five parameters to the de- 
sired accuracy of about 10 percent. One way to overcome 
this difficulty is to fit both the shear and the dielectric 
data with the same r c . This way was employed in our 
first two examples, PPE and DC704 (another diffusion 
pump oil). Then one obtains the other four parameters 
with reasonable accuracy from the shear data and e(0), 
Eoo and tjj from the dielectric data. One gets an ex- 
cellent fit, as one can see in Fig. 3. The parameters /3 
and /o = G Jg — 1 obtained turn out to be temperature- 
independent within their error bars. 



glass former 


PPE 


DC704 


TPE 


glycerol 


TPG 


T(K) 


254 


217 


250 


192.5 


192 


G (GPa) 


1.067 


1.602 


1.19 


3.49 


2.40 


■n (GPa s) 


0.0634 


0.216 


0.210 


9.44 


3.61 




2.37 


2.01 


2.2* 


2.2* 


7.0 


TD/r m 


8.1 


8.8 


5.9 


7.5 


19.0 


fo 


2.38 


2.46 


2.83 


4.61 


5.5 


$ shear 


0.64 


0.61 


0.53 


0.66 


0.64 




0.64 


0.61 


0.69 


0.667 


0.64 



TABLE I: Eshelby model parameters for five highly viscous 
liquids. The asterisk denotes set values. 

As it turns out, the value fo for Santovac 5P in Table 
I is a factor of 1.5 larger than the expectation of ir/2 of 
the Eshelby model. The deviation is clearly out of the 
error bar of about 10 %. This is in fact in agreement with 
shear compliance data in other molecular glass formers 19 , 
which show a similar tendency to higher values. Obvi- 
ously, the relaxational part of the recoverable compliance 
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FIG. 3: Model fits (continuous lines) of G"(lo) and e"(cu) 
(Roskilde data e ) for a diffusion pump oil at 254 K. 
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tends to be larger than the Eshelby expectation. We 
come back to this point in the discussion. 

The measurements of the Roskilde group£ comprise 
two other type-A glass formers (glass formers without 
a discernible Johari-Goldstein peak), TPE and DC704, 
which give similar results. However, in TPE it turned 
out to be necessary to allow for a sizable difference of the 
P in dielectrics and shear, setting the ratio r c /r m to the 
Eshelby expectation of 2.2 (see Table I). 

Glycerol turned again out to be an example where 
shear and dielectrics could be fitted with the same (3. 
Fig. 4 shows the fit of the shear data^. Their (3 agreed 
within the error bars with the one adjusted to dielectric 
datai^. Also, the times t c showed reasonable agree- 
ment, taking into account that these measurements come 
from different samples in different laboratories. 

As soon as one has a pronounced Johari-Goldstein 
peak, the relation r c = 2.2r m does not seem to hold any 
longer. In tripropylene glycol (Fig. 5), even the fit of the 
shear data alone already requires a considerably higher 
t c . This becomes even worse if one demands the same r c 
in dielectrics and shear. The fit in Fig. 5 was done with 
the same (3 for dielectrics and shear, but with r c = 6.9r m . 

IV. DISCUSSION AND CONCLUSIONS 

The comparison to experiment in the preceding section 
has raised more questions than it has answered. Insofar, 
the Eshelby model presented here is at present not more 
than one of the many alternatives in the field. It has, 
however, one distinct advantage, namely a direct con- 
nection to the elementary flow process. 

To illustrate this, let us consider the example of Santo- 
vac 5P or polyphenylene of Fig. 3. The molecule consists 
of five phenyl rings bonded by oxygen atoms. Such an 
oxygen bond is rather flexible. Therefore one should as- 
sess the typical molecular volume to one of these phenyl 
rings, as far as packing considerations are concerned. The 
density is about 1250 kg/m 3 , so the volume of one phenyl 




log(to/s" 1 ) 

FIG. 4: Fits of the shear data of Schroter and Donth 18 of 
glycerol. 



ring is 0.12 nm 3 . 

With this volume and the modulus of 1.06 GPa, the 
formation energy E a of the four-molecule strain defect of 
Fig. 1 according to eq. |T]) is 0.173 eV, corresponding 
to a temperature of 2004 K. The barrier height should 
be substantially lower because of the anharmonicity of 
the potential. This barrier height is too low; the barriers 
which cause the flow are centered around ksT g lnr m /ro, 
where To is the microscopic time of about 10 -13 s, which 
means a barrier height of 0.7 eV. On the other hand, the 
formation energy of the four-atom defect is high. In fact, 
it is too high; eq. (fTTj) tells us that to get a relaxational 
compliance of ir/2G (that is what we need according to 
eq. (|14|) requires a number of symmetric-potential four- 
molecule defects no per molecule given by 

jtt_ = n Q v 2 el 

2G 16fc B T' v ; 

Here v is four times the molecular volume v mo i and 
is 7r/3. Thus we find that we need about 0.05 four-atom 
defects per molecule to get the sample to flow. However, 
with a formation energy eight times higher than the glass 
temperature T g , it is more likely we will get less than 10 -3 
such defects per molecule. In fact, the four-atom defect of 
Fig. 1 seems a better candidate for the tunneling states 
which determine the low-temperature properties than for 
the flow defects; it even has the right order of magnitude 
of coupling to the strain. 

Looking at eq. (|2"3"|) , one realizes that the way out of 
this dilemma is to increase the volume of the core region. 
Let us consider a core region of hundred phenyl rings 
with a diameter of 2.8 nm with the same formation en- 
ergy as the four- atom defect. Then its spontaneous shear 
deformation e.; is only one fifth of 7r/3, a shear angle of 
12 degrees, but its coupling constant to the compliance 
is a factor of 25 higher. Therefore one needs a factor of 
25 less defects to make the sample flow, approaching the 
numbers which one expects on the basis of their Boltz- 
mann factor. 
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FIG. 5: Fits of G" (lo) and e"(u>) (Roskilde data 6 ) for tripropy- 
lene glycol, an example with a pronounced Johari-Goldstein 
peak. 
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A second important point is the high barrier. This 
cannot be a defect which reaches its saddle-point by 
only a deformation by 6j/2; there must be a more com- 
plicated rearrangement requiring an energy about four 
times higher. If the saddle point is indeed only reached 
by a high-energy rearrangement of the core region, one 
begins to understand why the relaxational compliance is 
higher than the prediction of eq. (114|) . There might be 
components of the motion which reverse after a while 
in following rearrangements, making the contribution of 
the defect to the flow smaller than expected. Imagine, for 
instance, that part of the rearrangement is a reorienta- 
tion of an interstitialiS. This would not contribute to the 
flow, but it would enhance the relaxational compliance. 
Another possibility is a molecular configuration change, 
which would also contribute to the compliance, but not 
to the flow, because it reverses after a while. The Johari- 
Goldstein peak is probably of this nature. But why the 
introduction of such a peak leads to a violation of the 
Eshelby lifetime condition, eq. ([5]) , as shown in the TPG 
example of Fig. 5, is not immediately clear. 

Another unsolved riddle (an old riddle^) is the short 
Debye relaxation time. The ratio of td and r m according 
to eq. (JTSJ) is 



T D _ ^Gr 3 H 



k B T 



(24) 



If we take ru from the molecular volume by 



Vmoi , this ratio for the phenyl ring in Santovac 5P at 254 
K should be 110, a factor of 14 higher than the one found 
in experiment. In terms of the radius m, this means that 
the radius must be only 40 percent of the one calculated 
from the molecular volume. One knows from NMR gradi- 
ent measurements of the molecular diffusion^! that the 
hydrodynamic radius is a bit smaller than the one cal- 
culated from the molecular volume, but not that much. 
Maybe the process which we characterize by Tr> has an- 
other and faster mechanism than the one considered by 
Debye, which might account for the fact that it is ac- 
companied by the relaxation of the energy and by the 
relaxation of the structu re 17 ' 22 . 

A third unanswered question is whether the structural 
rearrangements fully relax the molecular orientation, as 
assumed for simplicity in the derivation of equs. (j2l"j) 
and (f2"2")l . Their number is limited by the cutoff at V c , 
so it is conceivable that there remains a rest of dielec- 
tric polarization, which should appear as an additional 
Debye relaxation at the Debye relaxation time. In fact, 
dielectric measurements in the monoalcohols 2 ^ look as if 
most of the reorientation happened in one big Debye-like 
process long after the Maxwell time. 

We conclude that the Eshelby model for the highly 
viscous flow, though it raises more questions than it an- 
swers, provides a new way to tackle an old and difficult 
problem on a quantitative level. 
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